Diffusive-Ballistic Crossover and the Persistent Spin Helix 
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Conventional transport theory focuses on either the diffusive or ballistic regimes and neglects the 
crossover region between the two. In the presence of spin-orbit coupling, the transport equations 
are known only in the diffusive regime, where the spin precession angle is small. In this paper, 
we develop a semiclassical theory of transport valid throughout the diffusive - ballistic crossover of 
a special 5*17(2) symmetric spin-orbit coupled system. The theory is also valid in the physically 
interesting regime where the spin precession angle is large. We obtain exact expressions for the 
density and spin structure factors in both 2 and 3 dimensional samples with spin-orbit coupling. 

PACS numbers: 72.25.-b, 72.10.-d, 72.15. Gd 
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The physics of systems with spin-orbit coupling has 
generated great interest from both academic and practi- 
cal perspectives [l[ . Spin-orbit coupling allows forpurely 
electric manipulation of the electron spin [2,3.4.5.6], and 
could be of practical use in areas from spintronics to 
quantum computing. Theoretically, spin-orbit coupling 
is essential to the proposal of interesting effects and new 
phases of matter such as the intrinsic and quantum spin 
Hall effect 0JMMl3Il3II3[- 

Whilc the diffusive transport theory for a s yste m with 
spin-orbit coupling has recently been derived |l3llT3 ] , the 
analysis of diffusive-ballistic transport - where the spin 
precession angle during a mean free path is compara- 
ble to (or larger than) 2-7T - has so far remained con- 
fined to numerical methods jl5[ . This situation is exper- 
imentally relevant since the momentum relaxation time 
t in high-mobility GaAs or other semiconductors can 
be made large enough to render the precession angle 
4> — akpT > 2ir, where a, kF are the spin-orbit coupling 
strength and Fermi momentum respectively. The math- 
ematical difficulty in obtaining the crossover transport 
physics rests in the fact that one has to sum an infinite 
series of diagrams which, due to the spin-orbit coupling, 
are not diagonal in spin-space. In this paper we obtain 
the explicit transport equations for a the series of models 
with spin-orbit coupling where a special SU(2) symmetry 
has recently been discovered [l6j |. 

We first consider a two-dimensional electron gas with- 
out inversion symmetry for which the most general form 
of linear spin-orbit coupling includes both Rashba and 
Dresselhaus contributions: 
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where k x>y is the electron momentum along the [100] and 
[010] directions respectively, a, and (3 are the strengths of 
the Rashba, and Dresselhauss spin-orbit couplings and m 
is the effective electron mass. At the point a = (3, which 
may be experimentally accessible through tuning of the 



Rashba coupling via externally applied electric fields [|J , 
a new SU (2) finite wave- vector symmetry was theoreti- 
cally discovered [16|. The Dresselhauss [110] model, de- 
scribing quantum wells grown along the [110] direction, 
exhibits the above symmetry without tuning to a partic- 
ular point in the spin-orbit coupling space. At the sym- 
metry point, the spin relaxation time becomes infinite 
giving rise to a Persistent Spin Helix. The energy bands 
in Eq. [1] at the a = /3 point have an important shifting 
property: ej,(fc) = £f(fc + Q), where Q+ = Ama,Q- = 
for the 7i[R e D] model and Q x = Ama, Q y = for the 
Hm n] model. The exact 577(2) symmetry discovered in 
jl6| is generated by the spin operators (written here in a 
transformed basis as): 
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with Ckt,± being the annihilation operators of spin-up and 
down particles. These operators obey the commutation 
relations for angular momentum, [Sq,Sq] — ±25q and 
[Sq, Sq] = Sq . Early spin-grating experiments on GaAs 
exhibit phenomena consistent with the existence of such 
a symmetry point pV7| . 

In [l6| the spin-charge transport equations for the 
Hamiltonian Eq.pQ have been obtained in the diffusive 
limit in which akpT << 1. However the regions akpT ~ 
1 and akpT >> 1 are also experimentally accessible, and 
no theory is yet available to deal with these regimes. We 
now present the exact spin and charge structure factors 
at the exact symmetry point for any value of the param- 
eter akpT. 

We first obtain the spin and charge structure factors 
in the absence of spin-orbit coupling, but valid in both 
the r — > and in r — * oo regimes. One should think of 
the structure factor obtained this way as a generalization 
of the classic Lienhard formulas in the presence of disor- 
der. We then use a non-abelian gauge transformation 



2 



X 



-1 



FIG. 1: The sketch of the branch cut and the integral contour 
in the calculation of S(t,q). 



introduced in [16[ to obtain the structure factors for the 
spin-orbit coupling problem described above. 

We start by formulating the problem in the language 

33]. Assuming isotropic 



of the Keyldish formalism 
scattering with momentum lifetime r, the retarded and 
advanced Green's functions are: 



G"< A (k,e) = (e-H±-) 



(3) 



We introduce a momentum, energy, and position depen- 
dent charge-spin density which is a 2 x 2 matrix g(k, r, t). 
Summing over momentum: 



P(r,t) 



d 2 k 



g{k,r, t), 
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gives the real-space spin-charge density p(r, t) — n(r, t) + 
S l (r,t)(Ti, where n(r,t) and S % (r,t) are the charge and 
spin density and v — to/2-7t is the density of states in 
two-dimensions. p(r, t) and g(k 1 r,t) satisfy a Boltzman- 
type equation 

f4{f'£}+^ = -?+7< G ^>- < 5 > 

that we now solve for a free electron gas Hamiltonian. 
To obtain the spin-charge transport equations, we follow 
the general sequence of technical manipulations: time- 
Fourier transform the above equation; find a general so- 
lution for g{k,r,t) involving p{r,t) and the /c-dependent 
spin-orbit coupling; perform a gradient expansion of that 
solution (assuming d r << kp where kp is the Fermi 
wavevector) to second order; and, finally, integrate over 
the momentum. The formalism is valid even through the 
diffusive-ballistic boundary. For the diffusive limit, when 
r is small, we need to keep only the second order term in 
the gradient expansion which gives rise to the usual spin 
and charge propagator (iuj — Dq 2 )^ 1 . As r increases, we 
need to keep higher order terms in the gradient expan- 
sion to accurately describe the transport physics. The 
ballistic limit requires infinite summation over the gra- 
dient expansion. This can be easiest seen in the regime 



of zero spin-orbit coupling, in which the sums can be ex- 
actly performed. It is then fortuitous that our spin-orbit 
coupled problem can be mapped into a free electron plus 
disorder problem where we can obtain the structure fac- 
tor exactly. By Fourrier transforming in time we obtain 
the following recursive equation: 



iu>p(r 
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were O = lu + i/r and the n-th order term reads: 

g n (k, r, t) - d ri ...d rn ((-^)...(-^)(l)-g ( k , V, t) 

( ? ) 

where go(k, r, t) contains a term which fixes the momen- 
tum at the Fermi surface: 
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Since the initial Hamiltonian and the transport equa- 
tions are rotationally invariant we can assume propaga- 
tion only on [100] and with the use of the identities: 
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we can integrate over the Fermi surface angles to obtain 
the structure factor pole: 



S(u,q) = 



1 



iuj - 



(11) 



The correct interpretation of our structure factor requires 
consistently picking a branch of the square-root function 
in the denominator. We pick the branch cut along the 
positive a;-axis. The pole in the structure factor repre- 
sents the characteristic frequencies of the system: 



^1,2 
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which in the diffusive and ballistic limits reduces to the 
well known expressions: 
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Wi,2 ~ ±VFq 
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where D = v f t/2. The presence of only one (exponen- 
tially decaying) solution in the diffusive limit follows di- 
rectly from correctly treating the branch-cut singularity 
in our structure factor. It can then be seen that the ex- 
ponentially divergent solution ui « iDq 2 is a false pole of 
EqUU. 
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Although not of immediate interest to the present pa- 
per, we also present the structure factor for a bulk Fermi 
gas in the presence of disorder. With the density of states 

<2m) 3/2 E 1/2 

defined as v = - — ^ F the transport equation be- 
comes: 

f f f d<t> sin 6d6k 2 dk^^ . 
J J {2 *)*vt % 9n (W) 

I 



where g n and £1 are as before and n = lu+i/t. Rotational 
invariance allows us to take fcj = k z and we obtain: 



mk F ^ fv F q\ f 1 mkp Q ( 1 + l^n] 



Introducing the three-dimensional density of states at the 
Fermi surface, as well as a <5-function source term, the 
structure factor reads: 
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To see the diffusive pole we need to carefully expand the 
logarithm: 



0: p 



-^q 2 



(17) 



Which is the right diffusive pole in 3D. For the ballistic 
pole we solve the equation (the one below is valid for any 
t): 




FIG. 2: (a) The imaginary part and (b) the real part of S(t, q). 
We set t — 1. For both figures, from bottom to top, the curves 
are corresponding to a = 2.2, 2.6, 3, 3.4, 3.8, 4.2. 
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In the ballistic limit r — > oo the exponentials in the frac- 
tion are oscillating wildly and must be regularized. De- 
pending on on the regularization q — ► q + 0^ the charac- 
teristic frequencies are: 



u> = ±v F q 



(18) off-diagonal operators, such as S (x) — ipUx)ipi(x) 
and S + (x) = ip]-(x)ipi(x) are transformed: S~(x) — > 

exp(-iQ • r)S~(x), S+(x) -> exp(iQ • r)S+(x). Here 
Q is the shifting wavevector of the spin-orbit coupled 
Hamiltonian. Since in the gauge transformed basis, all 
three components of the spin and charge have the struc- 
ture factor derived above, in the original (experimentally 
measurable) basis, the S x and S y have the following form: 



(19) 



which are the ballistic poles. 

Having solved the free-Fermi gas case, we now add 
spin-orbit coupling at the special SU (2) symmetric point 
of the Persistent Spin Helix. Following [Ty] , we express 
the spin-orbit coupling Hamiltonian Eq.Q] in the form 
of a background non-abelian gauge potential HroD — 

tt~~ + -J—(k+ — 2maa z ) 2 + const, where the field strength 

2m 2m v ^ * > ° 

vanishes identically for a = f3. Therefore, we can elimi- 
nate the vector potential by a non-abelian gauge trans- 
formation: 4 r -[-(a; +) xJ) — * exp(i2max+)^i(x+,X-), 
^l(x+, X-) — ► exp(— i2max+) 1 § 1 i(x+, X-). Under this 
transformation, the spin-orbit coupled Hamiltonian is 
mapped to that of the free Fermi gas, but, while diagonal 
operators such as the charge n and S z remain unchanged, 
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iw — - + - 



(20) 



The above result represents the exact form factor for 
a spin-orbit coupled system valid everywhere from the 
diffusive to ballistic regimes. The Persistent Spin He- 
lix is clearly maintained for any values of T,a,Vf since 
S(u>, Q) — X/iu) which renders the spin life-time infinite. 

The transient grating experiments [l7.19] measure 
the W Fourrier transform of S(w,q), i.e. S(t,q) = 

J dte~ ZUJt S (u> , q) . S(ui,q) is analytic in the upper half 
complex plane. Thus, S(t, q) is zero for t < 0. For t > 0, 
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FIG. 3: The oscillation frequency Q in the imaginary part of 
S(t, q) as a function of a — vf \q ± Q\r. 



by selecting the integral contour as shown in fig.JT]), we 
obtain its real part and imaginary part as follows: 



Im(S(t,q)) a 



1 + a 2 



P 



2 Vx 2 - a 2 cos{^) 
a 7T x(x 2 — 1 — a 2 ) 



Re(S(t,q)) _ a 2 + cos{Vl + a 2 ±) 



1 + a 2 



(21) 



where a = vp\q±Q\r and P indicates the principal value 
of the integral. 

In Fig. (21) , we plot the real and imaginary part of 
S(t, q) for different values of a. In the figure, we set 



t = 1 and from bottom to top, the curves are correspond- 
ing to a = 2.2, 2.6, 3, 3.4, 3.8, 4.2. Although the real part 
is clearly an oscillating function of t with an oscillation 
frequency, ^ 1+a2 , the oscillation is not easily seen in the 
figure. However, the imaginary part has a much larger 
oscillation amplitude than the real part and the oscilla- 
tion becomes clear as increasing a, reflecting the ballistic 
nature of the sample. The oscillation frequency fl in the 
imaginary part is linearly dependent on a as shown in 
Fig.©. 

In this paper we have obtained the exact transport 
equations valid in the diffusive, ballistic, and crossover 
regimes of a special type of spin-orbit coupled system 
which enjoys an SU(2) gauge symmetry. We obtained 
the exact form of the structure factors, and found the 
dependence of the spin-density as would be observed in a 
transient-grating experiment. It would be interesting to 
work out the transport equations in the diffusive-ballistic 
regime in perturbation theory away from the Persistent 
Spin Helix. 
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